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Abstract 

We demonstrate the Casimir interaction between two ferromagnetic boundary surfaces using the 
dynamic atomic force microscope. The experimental data are found to be in excellent agreement 
with the predictions of the Lifshitz theory for magnetic boundary surfaces combined with the 
plasma model approach. It is shown that for magnetic materials the role of hypothetical patch 
potentials is opposite to that required for reconciliation of the data with the Drude model. 
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The Casimir effect whicli is a fluctuation-induced quantum phenomenon, is of much 
interest due to its promising multidisciplinary applications in nanotechnology, condensed 
matter physics, physics of elementary particles, and in gravitation and cosmology 
Many experiments on measuring the Casimir force between boundary surfaces made of 
different materials separated by a vacuum gap or a liquid have been performed in the last 
15 years 4t6|. It was shown that the magnitude of the Casimir force can be controlled by 

71, l8| , phase transitions 



using different boundary materials 



14 



17|. 



surfaces structured with nanoscale corrugations 

It is common knowledge that an unified description of both the van der Waals and Casimir 
forces is given by the Lifshitz theory [18] which characterizes the material properties by 
means of the frequency-dependent dielectric permittivity s{uj) and magnetic permeability 
Atfcu). The role of magnetic materials in the Casimir force has been studied theoretically 
19|-[30|. The interest stems from the possibility to obtain a repulsive Casimir force for 



application in micromachines. Using real magnetic materials 2l|, |25| did not validate the 
early results which used constant e and fi. As fi{i^) can be large only at ,^ < 10^ — 10^ Hz its 

m. 



entire contribution to the Lifshitz formula is through the zero Matsubara frequency 2^, 
For metals, the zero-frequency term is strongly influenced by the inclusion (Drude model 
approach) or neglect (plasma model approach) of the relaxation properties of free electrons 
^. Thus using fi provides another parameter to study the role of the relaxation properties 
of free electrons in the Casimir effect. Some experiments demonstrate strong disagreement 
between the measured data and theoretical predictions when the relaxation properties of 



electrons are taken into account for metals 131 



dielectrics 



43 q or the dc conductivity is included for 
13l . I36l439l|. The same data are found to be consistent with theory when 



the relaxation properties are neglected for metals or the dc conductivity of dielectrics is 
disregarded. Two other experiments |40. 



approach (see critical discussion in [42 



4l| are claimed to be in favor of the Drude model 
45l|). It was also hypothesized [461] that the effect 

341 ] in agreement with the 



of large patches might bring the experimental data of Refs. 
predictions of the Drude model approach (see also discussion in Ref. [35l] ) . 

In this Letter we describe demonstration of the Casimir force between surfaces of a plate 
and a sphere, both coated with ferromagnetic metal Ni, performed by means of dynamic 
atomic force microscope (AFM) using the frequency shift technique. The Lifshitz theory was 
generalized for the case of magnetic bodies in Ref. 19|], but till now was not unequivocally 
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verified experimentally. Note that measurements of the Casimir interaction between an Au- 
coated sphere and a Ni-coated plate confirmed the influence of magnetic properties on 
the Casimir force under an assumption that the plasma model approach is adequate (note 
that for a nonmagnetic metal interacting with a magnetic one the Drude model approach 
is not sensitive to magnetic properties and, for Au interacting with Ni, leads to almost 
the same results as the plasma model approach at the separations considered |4^). The 
advantage of Ni-Ni test bodies used here is that the magnetic properties significantly affect 
the Casimir force when both the plasma and Drude model approaches are used leading to 
considerably different results Using this property, we have unequivocally confirmed that 
the magnetic properties influence the Casimir force in accordance with predictions of the 
Lifshitz theory. The agreement is excellent with the plasma model approach, and the Drude 
model approach is excluded by our data at a 95% confidence level. We have also excluded 
any possible role of patch effects on the conclusions obtained. This opens opportunities 
for far-ranging applications of the magnetic Casimir effect in nanotechnology including the 



26-28 



realization of the Casimir repulsion through a vacuum gap 

Here we have used the same apparatus and cantilever preparation as in Refs. [35|, l47j . 
The gradient of the Casimir force was measured acting between a Ni-coated hollow glass 
microsphere of R = 61.71 ± 0.09 fim radius attached to the tip of a rectangular Si cantilever 
and a Si plate also coated with Ni. The thicknesses of Ni coating were 210 ± 1 nm and 
250 ± 1 nm on a sphere and a plate, respectively. The hollow sphere leads to higher resonant 
frequencies and mechanical Q-factors offering higher sensitivities. To promote adhesion of 
the Ni coating a 10 nm layer of Cr followed by 40 nm layer of Al was done first. The coatings 
were performed at 10"^ Torr. To achieve uniformity of Ni layers, the sample was rotated 
during evaporation of the metals. A coating rate ?5i 3 A/s was used. Both test bodies were 
cleaned following the special multi-step procedure to remove any attached adsorbates and 
debris (see jsSj for details). The cantilever was clamped in a specially fabricated holder and 
placed inside the vacuum chamber that was capable of reaching a pressure of 10"^ Torr by 
using mechanical, turbo and ion pumps. The Ni-coated plate was fixed on the top of the 
piezo with double sided vacuum adhesive tape. The movement of the piezo was calibrated 
by a fiber interferometer with 635.0 ± 0.3 nm laser source. The tube piezo contains a small 
magnet introduced by manufacturer which is not needed in this experiment. The initial 
magnetic field was measured to be ~ 100 Gs using a Hall probe gaussmeter. To prevent any 
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effects from this field, we inserted a piece of mu-metal between the top of the piezo tube 
and the Ni-coated plate. The residual magnetic field is below the detection resolution of 
0.1 Gs. Both the initial and residual fields do not depend on separations considered and do 
not contribute to the force gradient measured in our work. 

The d.„a,„ic .easu..e,„e.t .che.e the f.e,ue„e. .nodu.ation ™ode. a. ,„ Refs. HQ. 
was used. The directly measured quantity was the change of resonant frequency of the peri- 
odically driven cantilever which was detected by a phase locked loop system [48| (see details 
for our setup in Refs. 0, |47|). The driving frequency was kept near the resonance frequency 
of the cantilever to obtain the highest signal to noise ratio. The resonance frequency was 



detected with an optical interferometer j49|. To keep the interferometric cavity length be- 
tween the top of the cantilever and the end of the fiber fixed, we used a piezo above the 
cantilever, which was controlled by a proportional-integral-derivative feedback loop. This 
prevents errors in the sphere-plate separation distance a due to cantilever deflection from 
Casimir, F(a), and electrostatic, -Fei(a), forces. 

For small oscillations in the presence of an external force Ftot(a.) = Fc\{a) + F{a), the 
frequency shift Aw = uor — is expressed [ssl, |43] as Aw = —{ujQ/2k)Fl^^^{a). Here, Ur is the 
resonance frequency in the presence of Ftot, i^o is the natural resonance frequency, k is the 
spring constant of the cantilever, and a = -Zpiezo + zq (^^piczo is the plate movement due to the 
piezoelectric actuator and zq is the point of the closest approach between the two surfaces, 
which in our case is much larger than the separation on contact). The electric force between 
a sphere and a p , 
known function 

potential difference. In terms of the measured parameters, Au takes the form 



ate can be expressed as Fei(a) = X{a,R){Vi — Vq)^, where X{a,R) is the 
Vi are the voltages applied to the plate, and Vq is the residual 



Atu = -PiV,-Vof-CF'{a), (1) 

where C = Uo/{2k) and (3 = /3{zpiezo, Zo,C, R) = CX'{a,R), and allows determination of 
F'{a) from the measured Au. 

For electrostatic calibrations and measurements of Au, 11 different voltages in the range 
from -64.5 to 31.6 mV were applied to the Ni plate, while the sphere remained grounded. 
The plate was moved toward the sphere starting at the maximum separation of 2.3 /im 
and the corresponding Au was recorded at every 0.14 nm. Continuous triangular voltages 
at 0.01 Hz were applied to the tube piezo to move the plate toward the sphere. This set 
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calibration errors in ^piezo and dri 



C for different regions of a [35|, 



t 0.003 nm/s in the 
4^]. The parabohc 



of measurements was repeated three times. The small mechanical dri: 
sphere-plate separation distance was corrected as described in Refs. [sS . 
dependence of Aw on Vi was used to find Vq at each separation [isl. Note that Vq is separation 
independent indicating the lack of any adverse surface contaminants and the high quality of 
the measured data (see Fig. 1 where the best fit of Vq to the straight line leads to a slope 
equal to only 1.5 x 10~^mV/nm). The mean value of Vq = —17.7 ± 1 mV was found. The 
values of and C were found by a least fitting of the parabola curvature /3 in Eq. ([1]). The 
mean values are Zq = 221.1 ±0.4nm and C = 52.4 ± 0.16 kHzm/N (the errors are indicated 
at a 67% confidence level). Then we obtain /3 = {neQC R/ a^){l-2cia^ / - Ac2a^ / + . . .), 
where Ci and C2 are given in Ref. and eo is the permittivity of vacuum. The absence of 

t was confirmed by verifying the independence of zq and 
47|. After the values of zq and C were found, the measured 
Aw was converted into F!^o^{a) and the absolute separation distances were determined. 

Now the 33 values of the gradient of the Casimir force F'{a) at each a can be obtained from 
Eq. ([T]) by subtracting the contribution of -Fe'il*^)- They are shown in Fig. 2 at separations 
from 223 to 320 nm with a step of 2 nm. The statistical properties of the experimental data 
are characterized by the histogram shown in an inset to Fig. 2 at a = 250 nm. The histogram 
is described by Gaussian distribution with the standard deviation ap' = 0.92 /iN/m. The 
respective mean is F' = 74.17 fiN/m. The mean values of F'{a) as function of a (with a step 
of Inm) are shown as crosses in Fig. 3(a-d) where the arms of the crosses indicate the total 
experimental errors found at a 67% confidence level. The total errors are mostly determined 
by the systematic error which decreases from 1.2 to 1.0 /iN/m when a increases from 223 to 
250 nm and then to 0.8 /iN/m at a > 350 nm. The random error is equal to 0.18 /iN/m and 
does not depend on a. 



The experimental data for F'{a) were compared with predictions of the Lifshitz theory. 



The Lifshitz formula for magnetic materials 
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422l | was adapted for sphere-plate geometry 



using the proximity^ force approximation (this leads to < a/R, i.e., < 0.36% error at the 
shortest separation ^]) with the result 

^2 



^ ^ poo 

F'{a) = 2kBTR^ I qik±dk±Y^ 

1=0 



(2) 



Here, ks is the Boltzmann constant, T = 300 K is the temperature at the laboratory, 
qf = kj_ + ^f/c^, and = 2TikBTl/h with Z = 0, 1, 2, . . . are the Matsubara frequencies. 



The prime near the summation sign multiphes the term with / = by 1/2 and the sum 
with respect to a imphes a summation in the transverse electric [a = TE) and transverse 
magnetic {a = TM) polarizations of the electromagnetic field. The respective reflection 
coefficients are given by 

_ £iqi -h _ HiQi - k 

Siqi + h Hiqi + ki 

where kf = kl + eifxi^f/c^, and si = e{i^i), fii = ^ 

The dielectric permittivity ei was obtained from the tabulated optical data [51] for the 
complex index of refraction of Ni by means of the Kramers- Kronig relation. The data were 
extrapolated to zero frequency either by means of the Drude model or by means of the plasma 
model. The plasma frequency Up = 4.89 eV and the relaxation parameter 7 = 0.0436 eV 



have been used 
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52| . The magnetic properties of Ni were described by the static magnetic 



permeability /io = 110 (at T = 300 K for all / > 1 /i, = 1 [27]). 

The theoretical force gradients F'{a) were computed using Eqs. ([2]) and (j3l). The obtained 
values were corrected for the presence of surface roughness. The roughness profiles were 
investigated using an AFM and the r.m.s. roughness on the sphere and the plate was found 
to be 6s = 1.5 nm and 6p = 1.4nmLrespectively. At separations a > 223 nm this allows the 
use of the multiplicative approach 0, Q, ll?]. The theoretical results are shown in Figs. 2 
and 3(a-d) within different separation regions by the black and gray bands (their widths 
are deffned by the errors in the optical data) for the Drude and plasma model approaches, 
respectively. Note that at separations considered the difference between the predictions of 
the Drude and plasma models for F' is approximately proportional to a^^. As can be seen 
in Fig. 3, the Drude model approach is excluded by the data at a 67% confidence level over 
the region from 223 to 420 nm. The plasma model approach is in excellent agreement with 
the data over the entire measurement range. In Fig. 4 we plot the same measurement data 
for F', but with the total experimental errors determined at a 95% confidence level over 
the region from 223 to 350 nm (in an inset the interval from 300 to 350 nm is shown on an 
enlarged scale). As can be seen in Fig. 4, over the region from 223 to 350 nm the Drude model 
approach is excluded even at a higher, 95%, confidence level. At 420 nm < a < 1000 nm 
both the plasma and Drude model approaches are consistent with the data. It should be 
noted, however, that at large a the data become not informative. Thus, at a = 400, 550, 750, 
and 1000 nm the total relative experimental error is equal to 6%, 20%, 100%, and 321%, 
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respectively, whereas the respective differences between the predictions of the Drude and 
plasma model approaches are equal to 8.5%, 10%, 12%, and 14%. 

We emphasize that according to the Lifshitz theory the magnetic properties of Ni in 
Ni-Ni system significantly influence the gradient of the Casimir force in the framework of 
both theoretical approaches (they increase F' when the Drude model approach is used and 
decrease it if the plasma model approach is applied). Thus, our measurements unequivocally 
demonstrate the influence of magnetic properties on the Casimir force as is predicted by the 
Lifshitz theory combined with the plasma model approach. Of even greater importance 
is the fact that for two magnetic metals the Lifshitz theory predicts > where the 
Drude and plasma model approaches are indicated by the indices D and p (see Figs. 3 
and 4 where the black bands are above the gray). This is opposite to the case of two 
nonmagnetic metals where F^ < F^ 33N35|. From this it follows that the inclusion of effect 
of patches in the calculation for two magnetic test bodies is in principle incapable to bring 
the data in agreement with the Drude model approach because patches always lead to an 
additional attractive force. This finally proves that surface patches do not play any role 
in our experiments and confirms theoretical approach to patches suggested in Ref. [53[ and 



used in Ref. 



32|. 



In this experiment both interacting bodies are magnetic. Therefore it is necessary to 
analyze possible contribution of magnetic forces into the measurement results. This is done 
by considering two parallel Ni films of L^; x = 0.9^5 1.1 cm^ area and applying the 

5411. For films more than 150 nm 



general formulation of the proximity force approximation js, 
thickness the magnetization of each domain is perpendicular to the film surfaces, i.e., has 
only the 2;-component equal to iM^, where Mg = 435 emu/cm^ 
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57|. The magnetization 



of the first (1) and the second (2) films can be described by a function of two variables 
X, y). In order to obtain the pair of infinite films described by the periodic functions. 



(1,2), 



fl 2) 

we perform the periodic continuation of Mz ' (x, y) as odd function with the periods L 



and Ly and use the Fourier series 



oo oo 



klTX 



Mp)(x,,) = $:5:M(:,'^)sin^sin 



nny 



k=0 n=0 



(4) 



(1 2) 

Here, Mqq' = if the spontaneous magnetization is absent. 



Next, using the standard formalism developed in magnetic force microscopy 58|, |59|, one 



can calculate the magnetic field created by one Ni film and the magnetic force acting on the 



other film. Keeping in mind that the magnetic force between a pair of domains belonging to 
different films can be both attractive and repulsive, and that these domains have different 
size and are randomly arranged, the resulting magnetic force on a sphere is equal to zero 
under the condition that the spontaneous magnetization of at least one film is zero. 

To avoid the spontaneous magnetization of Ni films, we made them sufficiently thick and 
screened the weak environmental magnetic field in our experimental setup. If, however, 
there is some nonzero spontaneous magnetization of both films, the resulting gradient of the 
magnetic force acting on a sphere, although nonzero, is negligibly small. This is because 



the magnetic field near the center of a large film does not depend on z for z L^, Ly [60 1. 
For example, even for a fully magnetized film the gradient of the magnetic force acting on a 
sphere in the region of experimental separations is much less than 2 x 10~'^yuN/m, i.e., much 
less than the experimental error in the measurements of F'{a). 

To conclude, we have experimentally demonstrated that the magnetic properties of Ni 
influence the Casimir interaction as predicted by the Lifshitz theory combined with the 
plasma model approach. The Drude model approach in application to magnetic metals 
is excluded at a 95% confidence level. We have also shown that any hypothetical patch 
potential will only exacerbate the deviation from the Drude model approach. The obtained 
results allow realization of the Casimir repulsion through a vacuum gap which could lead to 
many potential applications in nanotechnology. 

Acknowledgments 

This work was supported by the DOE Grant No. DEF010204ER46131 (equipment, 
G.L.K., V.M.M., U.M.) and NSF Grant No. PHY0970161 (G.L.K., V.M.M., U.M.). G.L.K. 
and V.M.M. were also partially supported by the DFG grant BO 1112/21-1. 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. B 51, 793 (1948). 

[2] K. A. Milton, The Casimir Effect: Physical Manifestations of Zero-Point Energy (World 

Scientific, Singapore, 2001). 
[3] M. Bordag, G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Advances in the 

Casimir Effect (Oxford University Press, Oxford, 2009). 

8 



[4] G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Rev. Mod. Phys. 81, 1827 
(2009). 

[5] G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepanenko, Int. J. Mod. Phys. B 25, 171 

(2011). 

[6] A. W. Rodriguez, F. Capasso, and S. G. Johnson, Nature Photon. 5, 211 (2011). 
[7] F. Chen, G. L. KUmchitskaya, V. M. Mostepanenko and U. Mohideen, Phys. Rev. Lett. 97, 
170402 (2006). 

[8] S. de Man, K. Heeck, R. J. Wijngaarden, and D. lannuzzi, Phys. Rev. Lett. 103, 040402 
(2009). 

[9] G. Torricelh, P. J. van Zwol, O. Shpak, C. Binns, G. Palasantzas, B. J. Kooi, V. B. Svetovoy 
and M. Wuttig, Phys. Rev. A 82, 010101 (R) (2010). 
[10] G. Bimonte, Phys. Rev. A 78, 062101 (2008). 

[11] G. Bimonte, E. Cahoni, G. Esposito, L. Milano, and L. Rosa, Phys. Rev. Lett. 94, 180402 
(2005). 

[12] C.-C. Chang, A. A. Banishev, G. L. Khmchitskaya, V. M. Mostepanenko, and U. Mohideen, 

Phys. Rev. Lett. 107, 090403 (2011). 
[13] A. A. Banishev, C.-C. Chang, R. Castiho-Garza, G. L. KUmchitskaya, V. M. Mostepanenko, 

and U. Mohideen, Phys. Rev. B 85, 045436 (2012). 
[14] H. B. Chan, Y. Bao, J. Zou, R. A. Cirelh, F. Klemens, W. M. Mansfield and C. S. Pai, Phys. 

Rev. Lett. 101, 030401 (2008). 
[15] Y. Bao, R. Guerout, J. Lussange, A. Lambrecht, R. A. Cirelli, F. Klemens, W. M. Mansfield, 

C. S. Pai and H. B. Chan, Phys. Rev. Lett. 105, 250402 (2010). 
[16] H.-C. Chiu, G. L. Klimchitskaya, V. N. Marachevsky, V. M. Mostepanenko, and U. Mohideen, 

Phys. Rev. B 80, 121402(R) (2009). 
[17] H.-C. Chiu, G. L. Klimchitskaya, V. N. Marachevsky, V. M. Mostepanenko, and U. Mohideen, 

Phys. Rev. B 81, 115417 (2010). 
[18] E. M. Lifshitz, Zh. Eksp. Teor. Fiz. 29, 94 (1955) [Sov. Phys. JETP 2, 73 (1956)]. 
[19] P. Richmond and B. W. Ninham, J. Phys. C: Solid St. Phys. 4, 1988 (1971). 
[20] S. Y. Buhmann, D.-G. Welsch, and T. Kampf, Phys. Rev. A 72, 032112 (2005). 
[21] M. S. Tomas, Phys. Lett. A 342, 381 (2005). 

[22] S. J. Rahi, T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. Rev. D 80, 085021 (2009). 

9 



[23] Yu. S. Barash and V. L. Ginzburg, Usp. Fiz. Nauk 116, 5 (1975) [Sov. Phys. Usp. 18, 305 
(1975)]. 

[24] O. Kenneth, 1. Klich, A. Mann, and M. Revzen, Phys. Rev. Lett. 89, 033001 (2002). 
[25] D. lannuzzi and F. Capasso, Phys. Rev. Lett. 91, 029101 (2003). 
[26] F. S. S. Rosa, D. A. R. Dalvit, and P. W. Milonni, Phys. Rev. A 78, 032117 (2008). 
[27] B. Geyer, G. L. Khmchitskaya, and V. M. Mostepanenko, Phys. Rev. B 81, 104101 (2010). 
[28] G. L. Khmchitskaya, B. Geyer, and V. M. Mostepanenko, Int. J. Mod. Phys. A 25, 2293 
(2010). 

[29] N. Inui, Phys. Rev. A 84, 052505 (2011). 
[30] N. Inui, Phys. Rev. A 86, 022520 (2012). 

[31] R. S. Decca, E. Fischbach, G. L. Khmchitskaya, D. E. Krause, D. Lopez, and V. M. Mostepa- 
nenko, Phys. Rev. D 68, 116003 (2003). 

[32] R. S. Decca, D. Lopez, E. Fischbach, G. L. Khmchitskaya, D. E. Krause, and V. M. Mostepa- 
nenko, Ann. Phys. (N.Y.) 318, 37 (2005). 

[33] R. S. Decca, D. Lopez, E. Fischbach, G. L. Khmchitskaya, D. E. Krause, and V. M. Mostepa- 
nenko, Phys. Rev. D 75, 077101 (2007). 

[34] R. S. Decca, D. Lopez, E. Fischbach, G. L. Khmchitskaya, D. E. Krause, and V. M. Mostepa- 
nenko, Eur. Phys. J. C 51, 963 (2007). 

[35] C.-C. Chang, A. A. Banishev, R. Castillo-Garza, G. L. Khmchitskaya, V. M. Mostepanenko, 
and U. Mohideen, Phys. Rev. B 85, 165443 (2012). 

[36] F. Chen, G. L. Khmchitskaya, V. M. Mostepanenko, and U. Mohideen, Optics Express 15, 
4823 (2007). 

[37] F. Chen, G. L. Khmchitskaya, V. M. Mostepanenko, and U. Mohideen, Phys. Rev. B 76, 
035338 (2007). 

[38] J. M. Obrecht, R. J. Wild, M. Antezza, L. P. Pitaevskh, S. Stringari, and E. A. Cornell, Phys. 

Rev. Lett. 98, 063201 (2007). 
[39] G. L. Khmchitskaya and V. M. Mostepanenko, J. Phys. A: Math. Theor. 41, 312002 (2008). 
[40] A. O. Sushkov, W. J. Kim, D. A. R. Dalvit, and S. K. Lamoreaux, Nature Phys. 7, 230 (2011). 
[41] D. Garcia-Sanches, K. Y. Fong, H. Bhaskaran, S. Lamoreaux, and H. X. Tang, Phys. Rev. 

Lett. 109, 027202 (2012). 
[42] V. B. Bezerra, G. L. Khmchitskaya, U. Mohideen, V. M. Mostepanenko, and C. Romero, Phys. 

10 



Rev. B 83, 075417 (2011). 
[43] G. L. Klimchitskaya, M. Bordag, E. Fischbach, D. E. Krause, and V. M. Mostepanenko, Int. 

J. Mod. Phys. A 26, 3918 (2011). 
[44] G. L. Klimchitskaya, M. Bordag, and V. M. Mostepanenko, Int. J. Mod. Phys. A 27, 1260012 

(2012). 

[45] G. L. Khmchitskaya, M. Bordag, and V. M. Mostepanenko, larXiv:1208.T757V l. 
[46] R. O. Behunin, F. Intravaia, D. A. R. Dalvit, P. A. Maia Neto, and S. Reynaud, Phys. Rev. 
A 85, 012504 (2012). 

[47] A. A. Banishev, C.-C. Chang, G. L. Klimchitskaya, V. M. Mostepanenko, and U. Mohideen, 

Phys. Rev. B 85, 195422 (2012). 
[48] F. J. Giessibl, Rev. Mod. Phys. 75, 949 (2003). 

[49] B. C. Stipe, H. J. Mamin, T. D. Stowe, T. W. Kenny, and D. Rugar, Phys. Rev. Lett. 87, 
096801 (2001). 

[50] G. Bimonte, T. Emig, and M. Kardar, Appl. Phys. Lett. 100, 074110 (2012). 
[51] Handbook of Optical Constants of Solids, ed. E. D. Palik (Academic, New York, 1985). 
[52] M. A. Ordal, R. J. Bell, R. W. Alexander Jr., L. L. Long, and M. R. Querry, Appl. Opt. 24, 
4493 (1985). 

[53] C. C. Speake and C. Trenkel, Phys. Rev. Lett. 90, 160403 (2003). 
[54] B. V. Derjaguin, Kolloid. Z. 69, 155 (1934). 

[55] S. Hameed, P. Talagala, R. Naik, L. E. Wenger, V. M. Naik, and R. Proksch, Phys. Rev. B 
64, 184406 (2001). 

[56] W. H. Kraan and M. Th. Rekveldt, J. Magnetism & Magn. Mat. 5, 247 (1977). 

[57] O. V. Snigirev, K. E. Andreev, A. M. Tishin, S. A. Gudoshnikov, and J. Bohr, Phys. Rev. B 

55, 14429 (1997). 
[58] A. Wadas and P. Griitter, Phys. Rev. B 39, 12014 (1989). 

[59] D. Sarid, Scanning Force Microscopy. With Application to Electric, Magnetic and Atomic 
Forces (Oxford University Press, Oxford, 1994). 

[60] J. J. Saenz, N. Garcia, P. Griitter, E. Meyer, H. Heinzelmann, R. Wiesendanger, L. Rosen- 
thaler, H. R. Hidber, and H.-J. Giintherodt, J. Appl. Phys. 62, 4293 (1987). 



11 



300 400 500 500 700 800 900 1000 
a (nm) 



FIG. 1: The residual potential difference between a Ni-coated sphere and a Ni-coated plate as a 
function of separation. The mean value of Vq is shown by the gray line. 




FIG. 2: (Color online) Comparison between the nonaveraged experimental data for F' (gray dots) 
and theory (black and gray bands computed using the Drude and plasma model approaches, re- 
spectively). The inset shows the histogram for the measured F' at a = 250 nm. / is the fraction 
of 33 data points having the values of F' in the bin indicated by the respective vertical lines. The 
corresponding Gaussian distribution is shown by the dashed line. The black and gray vertical lines 
show the theoretical predictions of the Drude and plasma model approaches. 
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FIG. 3: Comparison between the experimental data for F' (crosses plotted at a 67% confidence 
level) and theory (black and gray bands computed using the Drude and plasma model approaches, 
respectively) . 



120 
100 
-5, 80 
^ 60 
40 



240 260 280 300 320 340 
a (nm) 

FIG. 4: Comparison between the experimental data for F' plotted at a 95% confidence level and 
theory (black and gray bands computed using the Drude and plasma model approaches, respec- 
tively) . 
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